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Abstract 
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MSG: 35L70, 35Q75, 53C25, 53C80, 53Z05. 

Keywords: Formal Linearization, Dispersionless PDEs, Integrability, Conformal 
Flatness, Einstein- Weyl geometry, Self-Duality. 



1 



1 Introduction 



Let 

F{X\U,U^^,U^^^J,...) =0 (1) 

be a partial differential equation (PDE), where u is a (scalar) function of the independent 
variables x^, . . . ,x"'. The formal linearization of (1) results upon setting u ^ u + ev, and 
keeping terms of the order e. This leads to a linear PDE for v, 

ipiv) = 0, (2) 

where ip is the operator of formal linearization, 

= Fu + Fu^^V^i + Fu^^^^jV^ixJ + . . . . 

(here V^i is the operator of total differentiation by x*, and etc are its iterations). 

Notice that, for nonlinear F, the formal linearization depends on the background solution 
u. For instance, the linearization of the dispersionless Kadomtsev-Petviashvili (dKP) 
equation, u^t — {uUx)x — Uyy = 0, reads as v^t — {uv)xx — Vyy = 0. The linearized equation 
(2) appears in a wide range of constructions and applications: 

• Stability analysis of a solution u is based on the investigation of the spectrum of the 
linearized equation (2): this goes back to Lyapunov. 

• Symmetries of the equation (1) correspond to solutions of the linearized equation 
of the form v = v{x\u,Uxi,UxixJ , ■ ■ ■), which are required to satisfy (2) identically 
modulo (1): this goes back to Lie. 

• Contact invariants of ordinary differential equations can be obtained from the 
Wilczynski invariants of linearized equations [13, 14]. 

• Generalized Laplace invariants of PDEs appearing in the context of Darboux in- 
tegrability can be obtained from the Laplace invariants of linearized equations 
[58, 2, 36, 38, 37]. 

• Integrability of ordinary differential equations can be seen from the structure of the 
differential Galois group of hnearized equations: it must be Abehan, see [44] and 
references therein. 

In general, coefficients of a linear PDE (in particular, its symbol) can be interpreted 
as tensor fields. The natural question arises: Can one read the integrability of a 
given PDE off the geometry of its formal linearization? In this paper we answer 
this question in the affirmative for the four particularly interesting classes of PDEs in 3D, 
namely 

Equations of type 1: 

{a{u))xx + {b{u))yy + {c{u))u + 2{p{u))xy + 2{q{u))xt + 2{r{u))yt = 0, (3) 
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this class was introduced in [16] in the context of the 'central quadric ansatz', see Remark 
at the end of Sect. 2. The corresponding integrability aspects were discussed in [28]. 
Equations of type 2: 



fllUxx + f22Uyy + hsUtt + '^fuUxy + '^fuUxt + 2/23%t = 0, 



(4) 



here the coefficients fij depend on the first order derivatives Ux,Uy^Ut only. Equations of 
this type can be called quasilinear wave equations, their integrability was analyzed in [7], 
see also references therein. 

Equations of type 3: 



Equations of this form are known as the dispersionless Hirota type, or Hessian type 
equations. Their integrability was studied in [47, 25]. 

Equations of type 4: 



here u is a two-component column vector, and A,B,C are 2x2 matrices. Equations of this 
form are known as systems of hydrodynamic type. Their integrability was investigated in 
[24]. It was demonstrated in [46] that 'generic' integrable equations of the types (4), (5) 
and (6) can be parametrized by generalized hypergeometric functions. 

Solutions to equations (3)-(6) carry a canonical conformal structure which can be 
defined as follows: the symbol of formal linearization is a symmetric (2,0)-tensor g'^ G 
r(S'^TM) on the base manifold M with coordinates {x^, x"^, x^) = {x, y, t), which depends 
on a finite jet of the background solution u (in case (6) we have to use the dispersion 
relation). This tensor is only defined up to multiplication by a non-zero factor, which 
makes our theory conformal. We will always assume to be non-degenerate, in this case 
the inverse (0,2)-tensor g G r{S^T*M) defines a metric^. In coordinates, g = gijdx^dx^ 
where gij is the inverse of the matrix of g'^. The conformal class [g] of the metric g is the 
key invariant responsible for the linearizability/integrability of the equations under study. 
Our main results can be summarized as follows: 

• Equations (3)-(6) are linearizable (by a transformation from the natural equivalence 
group) if and only if the corresponding conformal structures g are conformally fiat 
on any solution (an extra condition is needed in the case (6)). This provides a 
simple linearizability test based on the computation of the corresponding Cotton 
tensor (since all objects depend on a finite order jet of u, we do not rely upon 
smooth solvability of the equation (1), and carry out calculations formally using the 
geometric theory of PDEs). 

^ Since our approach to the dispersionless integrability is based on the method of hydrodynamic re- 
ductions [24], which generally reffers to hyperbolic systems, we assume the Lorentzian signature of g, but 
this requirement can be removed. 




(5) 




(6) 
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• Equations (3)-(6) are integrable by the method of hydrodynamic reductions [23] if 
and only if the corresponding conformal structures g are Einstein- Weyl on any solu- 
tion (again an extra condition is needed in the case (6)). Recall that the Einstein- 
Weyl structure consists of a symmetric connection D and a conformal structure g 
such that 

(a) the connection D preserves the conformal class: "^[g] = 0, 

(b) the trace-free part of the symmetrized Ricci tensor of D vanishes. 
In coordinates, this gives 

^k9ij = ^kQij, R{ij) = ^9ij, (7) 

where u = u^dx^ is a covector, R{ij) is the symmetrized Ricci tensor of D, and A is 
some function [10]. In fact one needs to specify g and uj only, then the first set of 
equations uniquely defines ro. We point out that for all examples considered in this 
paper the covector u is expressed in terms of g by the universal explicit formula 

Wfc = 2gkjV^s{g^^) + V^kilndeigij). (8) 

Note that in 3D this formula is invariant under the transformation g ^ Xg, u ^ 
u; + cilnA, which is characteristic of the Einstein- Weyl geometry. 

We recall that the Einstein- Weyl equations (7) are integrable by twistor-theoretic 
methods [33]. Thus, solutions of integrable PDEs carry 'integrable' geometry. Equiv- 
alently, one can say that second order dispersionless integrable systems in 3D (having 
non-degenerate symbol) can be viewed as reductions of the Einstein- Weyl conditions, 
which therefore play the role of a universal 'master-equation' [56]. Let us mention that 
relations of dispersionless integrable systems to the Einstein- Weyl geometry have been 
discussed in [55, 8, 9, 39, 18, 20, 21, 17, 31]. 

Given a class of integrable PDEs such as (4)-(6), the verification of the Einstein- Weyl 
conditions (7) can be a formidable task, primarily due to a rather intricate structure of the 
integrability conditions. A way to bypass computational difficulties is to use the result of 
Cartan [10] which says that the Einstein- Weyl property of a triple (©, g, u) is equivalent 
to the existence of a two-parameter family of surfaces which are null with respect to the 
conformal structure g (that is, tangential to the null cones of g), and totally geodesic in 
the Weyl connection D. In the context of dispersionless integrable systems, such surfaces 
are provided by the corresponding dispersionless Lax pairs: these consist of A-dependent 
vector fields X, Y which are required to commute modulo the equation, identically in the 
'spectral parameter' A (for some classes of PDEs discussed in this paper, e.g. (4) and (6), 
the existence of such Lax pairs was proved to be equivalent to the integrability by the 
method of hydrodynamic reductions [7, 24]). Note that X and Y may contain derivatives 
with respect to A. Taking integral surfaces of the distribution spanned by X, Y in the 
extended four-space with coordinates x, y, t, A, and projecting them down to the space of 
independent variables x, y, t, we obtain the required two-parameter family of null totally 
geodesic surfaces, see the proofs of Theorems 5, 9 for further details. Computationally, this 
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approach has an advantage, allowing one to avoid working with the full set of integrability 
conditions (the problem is to prove the existence of a Lax pair, and this is where the 
integrability conditions are needed). In this approach, the key object within the triple 
(D, g, uj) is the connection D, which is uniquely specified by the given two-parameter 
family of null totally geodesic surfaces. Both Einstein- Weyl conditions (7) will be satisfied 
automatically. 

Example. As an illustration let us consider the dKP equation, 

Uxt ~ {UUx)x ~ Uyy = 0. 

It was demonstrated in [18] that the corresponding Einstein- Weyl structure is provided 
by the conformal metric g = Adxdt — dy"^ + Audt'^ and the covector u = —Auxdt. The 
dispersionless Lax pair is given by 

X = dy- Xdx + u^dx, Y = dt- {X"^ + u)dx + {u^X + Uy)dx, 

one can verify that the condition [X, Y] = is equivalent to dKP. Projecting integral 
surfaces of the distribution spanned by X, Y to the space of independent variables x, y, t, 
one obtains a two-parameter family of surfaces which are null with respect to g, and 
totally geodesic in the Weyl connection D specified by g and u. 

Our main results relating linearizability /integrability to geometry of formal lineariza- 
tions are proved in Sect. 2-5 (we find it more convenient to treat the above four classes 
separately: explicit forms of the corresponding linearizability/integrability conditions are 
rather different). Known integrable equations of types (3)-(6) provide an abundance of 
Einstein- Weyl structures parametrized by elementary functions, elliptic functions, mod- 
ular forms and Painleve transcendents. Some further examples are collected in Sect. 6. 
Most of the examples of Einstein- Weyl structures exhibited in this paper are apparently 
new (otherwise a reference is given). 

In Sect. 7 we discuss geometric aspects of integrability of second order dispersionless 
PDEs in 4D, indicating that the associated conformal structures must be self-dual. This is 
in agreement with the fundamental fact that the Einstein- Weyl equations are reductions 
of the equations of self-duality [35, 8, 9]. 

In calculations of the Cotton tensor and the Einstein- Weyl conditions we used symbolic 
packages of Maple. We shall omit the unnecessary lengthy formulae from the text. All 
results of this paper are local, though we perceive that global versions may be available 
through the twistor theory. 

2 Equations of type 1 

In this section we consider equations of the form (3), 

ia{u))xx + {b{u))yy + {c{u))tt + 2{piu))xy + 2{q{u))xt + 2{r{u))yt = 0. 

Their integrability was investigated in [28] based on the method of hydrodynamic re- 
ductions. This approach requires the existence of multi-phase solutions of the form 
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u = u{R^, . . . , K^), where the phases R^{x,y,t) satisfy a pair of consistent systems of 
hydrodynamic type, 

R\ = X\R)Rl., R'y = ^\R)Rl., 

see [30, 23]. The requirement of the existence of 'sufficiently many' solutions of this type 
imposes strong constraints on the coefficients of the equation, and provides an efficient 
classification criterion. To formulate the classification result we introduce the symmetric 
matrix 

fa' p' q'\ 
V{u) = I p' b' r' , 
y q' r' c' I 

where prime denotes differentiation by u. The classification is performed modulo (com- 
plex) linear changes of the independent variables x, t, as well as transformations u 
ip{u), which constitute the equivalence group of our problem. 

Theorem 1 [28] Equation (3) is integrable by the method of hydrodynamic reductions if 
and only if the matrix V{u) satisfies the constraint 

V" = {IndetVyV + kV, (9) 

for some scalar function k. Modulo equivalence transformations, this leads to the five 
canonical forms of nonlinear integrable models: 

Uxx + Uyy - (ln(l - e''))yy " (ln(l " e''))tt = 0, 

Uxx + Uyy + (e")tt = 0, 
(e" - u)^^ + 2u^y + (e")H = 0, 

U'Xt ~ {UUx)x ~ Uyy = 0, 
{U'^)xy + Uyy + 2U^t = 0. 

Examples 2 and 4 o'^e the familiar Boyer-Finley (BE) and the dKP equations. 

We point out that the constraint (9), which implies V" G span{V.,V'}, means that 
the 'curve' V{u) lies in a two-dimensional linear subspace of the space of 3 x 3 symmetric 
matrices. The classification of normal forms of such linear subspaces, which is equivalent 
to the classification of pencils of conies, leads to the five canonical forms of Theorem 1. It 
was pointed out by D. Calderbank that equations of the form (3) are related to generalized 
Nahm equations with the gauge group SDif fiT?). In this language, the five canonical 
forms of Theorem 1 correspond to the five types of generalized Nahm equations obtained 
in [8]. For the dKP equation this correspondence was explicitly demonstrated in [16]. 

The linearized equation (3) is 

o!{u)v^^ + b'{u)vyy + c!{u)vtt + 2p\u)v,j.y + 2q'{u)v^t + 2r\u)vyt H = 0, 

where dots denote terms with lower order derivatives of v. Its symbol defines a conformal 
structure g = gij{u)dx^dx^ where {x^,x'^,x^) = {x,y,t), and the matrix gij is the inverse 
of V. Our first result is as follows. 
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Theorem 2 Equation (3) is linearizable by a transformation from the equivalence group 
if and only if the conformal structure g is conformally fiat on any solution. 

Proof: 

The condition responsible for conformal flatness in three dimensions is the vanishing 
of the Cotton tensor, 

Vr{Rpq - -^Rdpq) = q{Rpr " ^-Rfl'pr); (10) 

where Rpq is the Ricci tensor, R is the scalar curvature, and V denotes covariant dif- 
ferentiation in the Levi-Civita connection of g. Calculating (10) and using (3) and its 
differential consequences to eliminate all higher order partial derivatives of u containing 
differentiation by t more than once, we obtain expressions which have to vanish identically 
in the remaining higher order derivatives of u (without loss of generality we will assume 
that c{u) = u, this can be achieved by a transformation from the equivalence group). 
Requiring the vanishing of coefficients at the remaining derivatives of u, we obtain that 
all entries of the matrix V must be constant, which leads to linear equations. This flnishes 
the proof of Theorem 2. 

It turns out that conformal structures corresponding to all flve integrable models from 
Theorem 1 satisfy the Einstein- Weyl property. In fact, this follows from the construction 
of [8] which provides Einstein- Weyl structures from solutions of the gauge fleld equations 
with the gauge group SDif f{T?) modelled on Riccati spaces; our goal here is to present 
the explicit formulae. 

Equation 1: u^^ + Uyy - (ln(l - e'^))yy - (ln(l - e'^))tt = 0. 
Conformal structure^: g = dx^ + (1 — e^)dy'^ + (e^" — l)dt'^. 
Covector: uj = ^^^u^dx — Uydy + Utdt. 

Equation 2: + Uyy + {e^)tt = (BP equation). 
Conformal structure: g = dx"^ + dy"^ + e~^dt^. 
Covector: u = —u^dx — Uydy + Utdt. 

This Einstein- Weyl structure was obtained in [55], see also [39]. 

Equation 3: (e" - u)^^ + 2u.j,y + (e'')tt = 0. 
Conformal structure: g = 2dxdy + (1 — e^)dy'^ + e~^dt^. 
Covector: uj = —u^dx + {2e^Ux — Uy)dy + Utdt. 

Equation 4: Uxt — {uUx)x — Uyy = (dKP equation). 
Conformal structure: g = Adxdt — dy"^ + Audt^. 
Covector: uj = —Au^dt. 

This Einstein- Weyl structure was obtained in [18]. 

Equation 5: {v^)xy + Uyy + 2uxt = 0. 

^Here and in what follows we choose the 'simplest' representative metric g within the conformal class 
[g] . In all cases it is given by either the inverse or the cofactor matrix of . 
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Conformal structure: g = 2dxdt + dy^ — 2udydt + v^dt^. 
Covector: uj = 2uxdy + 2{uy — uUx)dt. 



In general, the corresponding conformal structures g and covectors oj can be written in 
terms of the matrix V{u) as follows: 



g = {dx dy dt)V-^ \ dy \ , to = 2{dx dy dt)V^W % - d{\ndeiV). (11) 





They satisfy the Einstein- Weyl equations (7) if and only if V{u) satisfies the integrability 
condition (9). Setting {x,y,t) = {x^,x'^,x^) one can represent the components of u = 
ojkdx^ by the formula (8). It turns out that exactly the same formula holds for all other 
classes of dispersionless PDEs discussed in this paper. The covector oj is related to the 
symbol of formal linearization via the identity 

where V* = g^^Vk, = g^^oJk, and V denotes covariant differentiation in the Levi- 
Civita connection of the metric gij. Note that the right hand side of this identity can be 
interpreted as a special case of the Weyl wave operator of weight zero [20]. The second 
main result of this section is 



Theorem 3 Equation (3) is integrable by the method of hydrodynamic reductions if and 
only if the corresponding conformal structure g is Einstein- Weyl on any solution, with the 
covector u given by (8). 

Proof: 

Given an equation of the form (3), the conformal structure of its formal lineariza- 
tion is g = [dx dy dt)V~^{dx dy dty. We will seek a covector u in the form u = 
{dx dy dt)T{ux Uy UtY where T{u) is an unknown 3x3 matrix depending on u. Imposing 
the Einstein- Weyl equations (7) and using (3) to eliminate the second order derivative Uu 
(no higher order derivatives of u will occur in this calculation), we obtain a set of relations 
which have to vanish identically in the remaining partial derivatives of u. Thus, equating 
to zero terms at the second order derivatives of u, we obtain T in terms of V, 

T = 2V~^V' - {IndetVyE, 

where E is the 3x3 identity matrix. This is equivalent to the formula (8) for u. The 
remaining terms vanish identically if and only if V satisfies the constraint (9). This finishes 
the proof of Theorem 3. 

Remark 1. It was assumed in the proof of Theorem 3 that u depends linearly on the first 
order derivatives of u. One can show that this assumption is unnecessary: even a rather 
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general requirement that u depends on some finite order jets of u is already sufficiently 
restrictive, and leads to the formula (8) for u. 

Remark 2. Equations from Theorem 1 possess implicit solutions u{x,y,t) of the form 

{x,y,t)M{u){x,y,tf = 1, (12) 

where M{u) is a 3 x 3 symmetric matrix of u. The level surfaces of such solutions, u = 
const, are central quadrics in the space of independent variables x, y, t. This construction 
is known as the central quadric ansatz [53, 16]. The equations for M{u) is 

M' = cMVM/VdetM, 

c = const. It was demonstrated in [53, 16] that in the cases of BF and dKP, this equation 
reduces to Painleve transcendents P3-P1. It was shown in [28] that other integrable 
models from Theorem 1 lead to the remaining Painleve equations P6-P4, with the full 
P6 corresponding to the first equation. Thus, we obtain a whole variety of Einstein- Weyl 
structures parametrized by Painleve transcendents. 

3 Type 2: quasilinear wave equations 

In this section we discuss geometric aspects of quasilinear wave equations (4), 

fllUxx + f22Uyy + /sS^M + ^fuU^y + ^fi^Uxt + '2f23Uyt = 0, 

here the coefficients fij are assumed to be functions of the first order derivatives u^, Uy, Ut 
only. PDEs of this type were investigated in [7] based on their correspondence with confor- 
mal structures in projective space. It was pointed out that the moduli space of integrable 
equations is 20- dimensional. In was shown in [46] that 'generic' integrable equations of 
the form (4) can be parametrized by generalized hypergeometric functions. We recall 
that the class of quasilinear wave equations is invariant under the group GL(4) of linear 
transformations of the variables x*, u, where (x^, x^, x^) = [x, y, t). These transformations 
constitute the natural equivalence group of the problem. 
The linearized equation has the form 

fij'^x'i-xi + ■ ■ ■ = 0, 

and defines a conformal structure g = gijdx'^dx^ where the matrix of gij is the inverse of 
fij. Here dots denote terms with lower order derivatives of v. Our first result is as follows. 

Theorem 4 Equation (4) is linearizable by a transformation from the equivalence group 
GL{4) if and only if the conformal structure g is conformally flat on any solution. 

Proof: 
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Let us first recall, following [7], the linearizability conditions for equations of the form 
(4). Since the structure of these conditions is the same in any dimension, we will consider 
the general n- dimensional case, 

where fij are functions of the first order derivatives u^k only, i, j, k = 1, . . . ,n. Setting = 
u^k we will write down a system of differential constraints for fij{p) which are necessary 
and sufficient for the linearizability of the equation under study by a transformation from 
the equivalence group GL{n + 1). Let us introduce the object 

where 

^ (n+2)(l-7i) {^Pkfij " ''T-dpjfik) , 
= inj){n-l) ((^ + ^)dpkfij - 2('^ + l)dpjik) ■ 

Then the linearizability is equivalent to the following two conditions: 
• o-ijk = 0, 

Geometrically, these conditions are equivalent to the existence of a flat connection V (in 
p-coordinates) with Christoffel symbols T^j/^ = Sj6l + SkSj such that Vkfij = Ckfij, see [7] 
for more details. 

Let us now require that the conformal structure g = gijdx^dx^ is conformally fiat on 
any solution. Calculating the Cotton tensor (10) (now we set n = 3), and using (4) and 
its differential consequences to eliminate all higher order partial derivatives of u which 
contain differentiation by t more than once, we obtain a complicated expression which has 
to vanish identically in the remaining higher order derivatives of u. In particular, requiring 
that coefficients at the remaining fourth order derivatives vanish identically (no higher 
order derivatives of u will occur in this calculation), we obtain the first linearizability 
condition, aijk = 0. 

There are two ways to proceed: collecting terms at the lower order derivatives of 
u one can obtain the second set of linearizability conditions (this, however, leads to 
quite complicated calculations). Another way is to point out that the condition aij^ = 
alone is already sufficiently restrictive [7, 51, 1], and implies that the PDE in question is 
either linearizable, or reducible to the equation for minimal hypersurfaces in a (pseudo) 
Euclidean space, 

[{Vuf - l]Au - {Vu)H{Vuy = 0, 

where Vm = (m^i, ...,Mx") is the gradient of u, A is the Laplacian, and H is the Hessian 
matrix of u. In the three-dimensional case we arrive at the equation 

{ul + Ui- l)Uxx + {ul + Uf- l)Uyy + {ul + Ul- l)Utt 

-2{UxUyUxy + UxUtUxt + UyUtUyt) = 0. 
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To complete the proof it remains to point out that the corresponding conformal structure 
g is not conformally flat on any solution, furthermore, the equation itself is not linearizable 
(in fact, not even integrable for n > 3). This finishes the proof of Theorem 4. 

Remark. This proof, and the proof of Theorem 2, generalize to any dimension n > 3, 
with the only change that one needs the Weyl tensor of conformal curvature instead of the 
Cotton tensor. For non-linearizable differential equations, the requirement of conformal 
flatness singles out a subclass of exact solutions of a given PDE, see [41]. 

Our next goal is to prove that conformal structures corresponding to formal lineariza- 
tions of integrable equations of the form (4) give rise to the Einstein- Weyl geometry. Let 
us begin with examples of known integrable PDEs. 

Example 1. The equation UxUyt + UyUxt + UtUxy = constitutes the Euler-Lagrange 
equation for the Lagrangian density u^UyUt which was obtained in [26] in the classification 
of first order integrable Lagrangians. 

Conformal structure: g = [u^dx + Uydy + UtdtY — 2uldx'^ — 2Uydy'^ — 2u^dt'^. 
Covector: u = -4^ dx - 4^ dy - 4^ dt. 

UyUt UtUx UxUy 

Example 2. The equation [uypiut))^ + {uxp{ut))y + {uxUyp'{ut))^ = constitutes the 
Euler-Lagrange equation for the Lagrangian density UxUyp{ut), which can be viewed as 
a deformation of Example 1. In this case the integrability conditions reduce to a single 
fourth order ODE for p, 

nil I 2 II o '2\ 2 w2 ^ I II III , o 1'^ III n '2 //2 ^ 

p [p p — 2pp ) — p p + 2pp p p +8p p — 9p p =0. 

It was shown in [27] that the general solution to this ODE is a modular form of weight 
one and level three known as the Eisenstein series -£^1,3. 
Conformal structure: 

g = {p'uxdx + p'uydy + pdtY — 2p''^uldx'^ — 2p''^u^dy'^ — 2p^dt^ — 2pp"uxUydxdy. 

Covector: 

= 2 (^^ p(Jp"-2p'^) ^ ~ ut^ dx + 2 (^^ p{pp"^^2p'^) ^ ^y^tt - 2^ m^) dy 

This structure is Einstein- Weyl if and only if p satisfies the above fourth order ODE. 

The second main result of this Section is as follows: 

Theorem 5 Equation (4) is integrable by the method of hydrodynamic reductions if and 
only if the corresponding conformal structure g is Einstein- Weyl on any solution, with the 
covector u given by (8). 
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Proof: 

We will give two proofs of this result. The first one is computational, based on the 
explicit calculation of the Einstein- Weyl constraints, and the integrability conditions as 
derived in [7]. The second proof utilises the fact that any PDE of the form (4) possesses 
a dispersionless Lax pair [7]. We demonstrate that the existence of a Lax pair implies 
that the Weyl connection 3, specified by the conformal structure g = Qijdx^dx^ and 
the covector (8), possesses a two-parameter family of null totally geodesic surfaces, the 
property known to be characteristic of the Einstein- Weyl geometry [10]. 

The first proof can be summarized as follows. Given an equation of the form (4), the 
conformal structure of its formal linearization is g = gijdx^dx\ where gij is the inverse of 
fij. We will seek a covector u = ojkdx^ in the form Uk = Tl^'u^i^j, where T^-' are certain 
functions of the first order derivatives of u. Imposing the Einstein- Weyl equations (7) 
and using (4) and its differential consequences to eliminate all higher order derivatives 
of u that contain differentiation by t more than once (maximum third order derivatives 
of u will occur in this calculation), we obtain a set of relations which have to vanish 
identically in the remaining partial derivatives of u. Thus, equating to zero terms at 
the third order derivatives of u, we obtain the expression (8) for u. The remaining terms 
vanish identically if and only if the coefficients fij satisfy the set of integrability conditions 
as derived in [7]. This finishes the proof of Theorem 5. 

Let us now give a somewhat more conceptual (as well as less computational) demon- 
stration that the integrability implies Einstein- Weyl property. It is based on the fact that 
any integrable equation of the form (4) possesses a dispersionless Lax pair of the form 

St = f{S,, '^Xl Uy, Itf), Sy g(ySx, Ux, Uy, Uf) . (13) 

This means that the consistency condition. Sty = Syt, is equivalent to the equation (4). 
Lax pairs of this form are known to arise as dispersionless limits of solitonic Lax pairs in 
2 + 1 dimensions [57]. Let us first outline the general construction which leads from the 
Lax pair (13) to totally geodesic null surfaces of the Weyl connection ©. Differentiating 
(13) by X and setting 5*^. = A, = a, Uy = b, Ut = c we obtain 

At = fxK + fattx + fbbx + fcCx, \ = 9\K + QaO-x + Qb^x + QcCx- (14) 

With this system we associate the vector fields 

d d Odd d 

X = — - /a— + {fattx + fhhx + fcCx)-Qy^, ^ ^ 'dy~ ^^dx ^ ^ ^^^"^ ^"^''''aA' 

which live in the extended four-dimensional space with coordinates x,y,t,X. Note that 
the compatibility condition Xty = Xyt is equivalent to the commutativity of these vector 
fields: [X, Y] = 0. The geometry behind this construction is as follows. Let us consider 
the cotangent bundle Z of the solution u{x, y, t), with local coordinates (x, y, t, Sx, Sy, St). 
Equations (13) specify a four- dimensional submanifold C Z parametrised by x, y, t and 
A. The compatibility of the equations (13) indicates that this submanifold is coisotropic. 
The vector fields X, Y are nothing but generators of the kernel of the restriction to of 
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the symplectic form dS^ /\dx + dSy Ady + dSt A dt. Equations (14) mean that the vectors 
X, Y are tangential to the hypersurface of defined by the equation A = A(x, y, t). 

Projecting the two-parameter family of integral surfaces of the distribution spanned 
by X, Y to the space of independent variables x, y, t we obtain a two-parameter family 
of null totally geodesic surfaces of the Weyl connection D. To see this we first project X 
and Y . This gives two vector fields 

^ _ d d ^ _ d d 
ot OX oy ox 

which commute if and only if A satisfies the equations (14). It remains to show that X 
and Y form a null distribution (that is, tangential to the null cones of g), and that the 
covariant derivatives 'S^j^X, Dj^y, DyX, D^y belong to the span of X,Y. Equivalently, 
one can introduce the covector 6 = dx + g\dy + f\dt which annihilates X, Y, and verify 
that 6 is null, and that D^^ A ^ = HyO A6 = 0. This follows from the equations satisfied 
by the functions /(A, a, b, c) and g{\, a, b, c) as derived in [7]: 

fa = 2kfi2 + gxkf 22 + fxikf23 " P), 9a = -2kfi3 - fxkfss - g\{kf 23 + p), 

fb = 9b = ~kf23 + P, (15) 

fc = kf23+P, 9c=-kfs3, 

where p{X, a, b, c) and k{X, a, b, c) are yet another two auxiliary functions. Furthermore, 
fx and gx satisfy the relation 

fii + f229l + /33/A + 2/i2^7a + 2/13/A + 2f23fx9x = 0, (16) 

which means that the covector 6 is null. To close the system (15) - (16) one proceeds as 
follows. Calculating the consistency conditions for Eqs. (15), fab = fba, 9ab = 9ba^ etc, six 
conditions altogether, and differentiating the relation (16) by a, 6, c and A, one obtains 
ten relations which can be solved for fxx,9xx and the first order derivatives of k and p. 
Modulo these relations, the conditions V>^6 A 6 = Dy^ A ^ = are satisfied identically. 
This finishes the proof of Theorem 5. 

4 Type 3: dispersionless Hirota equations 

In this section we discuss geometric aspects of PDEs (5) of the dispersionless Hirota type, 

which were investigated in [47, 25], revealing a remarkable correspondence with hypersur- 
faces of the Lagrangian Grassmanian. Geometric aspects of GL{2) structures associated 
with such equations were studied in [52]. In was shown in [46] that the 'generic' in- 
tegrable equation of the form (5) can be parametrized by generalized hypergeometric 
functions. Recall that equations of the form (5) are invariant under the group Sp{6) of 
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linear symplectic transformations of the variables x'^,u^i, where {x^,x'^,x^) = {x,y,t). 
These transformations constitute the natural equivalence group of the problem. 
The linearized equation has the form 

and defines a conformal structure g = Qijdx^dx^ where the matrix of Qij is the inverse of 
Fij = dF/du^i^j ■ Our first result is as follows. 

Theorem 6 Equation (5) is linearizable by a transformation from the equivalence group 
Sp{Q) if and only if the conformal structure g is conformally flat on any solution. 

Proof: 

Solving for Uu (this can be done without restricting the generality due to the presence of 
a large equivalence group), we can rewrite (5) in the form 

Utt = f{u )• (17) 

The corresponding linearized equation is 

"^tt fuxx'^xx "I" fuxy'^xy "I" fuyy'Vyy "I" fu^t'^xt "I" fuyt'^yti 

with the associated conformal structure (with upper indices) defined by the matrix 

2fuxy fuyy 2'^"yt 
V 2.^"^t 2f^yt ~1 / 

We require this structure to be conformally fiat for any background solution u{x,y,t). 
Calculating the Cotton tensor (10) and using (17) and its differential consequences to 
eliminate all higher order derivatives of u containing differentiation by t more than once, 
we obtain a complicated expression which has to vanish identically in the remaining higher 
order derivatives of u (maximum fifth order derivatives will appear in this calculation). In 
particular, requiring that coefficients at the remaining fifth order derivatives of u vanish 
identically, we obtain nine second order differential constraints for /: 

fuXxfuxtUxt ~l~ fuxxUxX 0) fuyyfuytUyt ~^ f UyyUyy 0, 

fuxtfuxtUxt '^fuxtUxX 0, fuytfuytUyt ~^ "^fuytUyy 0, 

fuytfuxtUxt ~^ '^ifuxtfuxtUyt ~^ fuxtUxy ~l~ fuytUxx) 0, 

fuxtfuytUyt ~l" '^ifuytfuxtUyt ~^ fuytUxy ~^ fuxfUyy) 0, (^S) 

fuxyfuxtUxt ~^ '^fuxxfuxtUyt ~^ "^fuxxUxy 
fuxyfuytUyt ~l" "^fuyyfuxtUyt ~^ fu yy'^xy 7 
fuyyfuxtUxt ~l~ fuxxfuytUyt ~^ fuxy fuxtUyt ~^ f Uxx^yy + fux yUxy 

14 



It was shown in [50, 6, 12] that these relations characterise symplectic Monge- Ampere 
equations, that is, PDEs (5) such that the left hand side F can be represented as a linear 
combination of all possible minors of the Hessian matrix of m, 

('^xx Uxy Uxt \ 
^xy ^yy ^yt I • 
Uxt Uyt Utt I 

Symplectic Monge- Ampere equations and differential constraints (18) have a clear geo- 
metric interpretation. Let us consider the Lagrangian Grassmannian A^ which can be 
(locally) parametrised by 3 x 3 symmetric matrices U . Minors of U define the Pliicker 
embedding of A^ into projective space P^^ . We will identify A^ with the image of this 
projective embedding. Symplectic Monge- Ampere equations can be viewed as hyperplane 
sections of A^ C P^^ . We point out that differential constraints (18) can be repre- 
sented in compact form as 

d^f = 2ao{dfdUxy - duxtduyt) + 2ai{dfduyy - {duytf) + 2a2{dfduxx - {dUxtY) 

(19) 

+2bo{duxxduyy - {duxy)"^) + 2bi{duxtduxy - duytduxx) + 2h2{duytduxy - duxtduyy) 

indeed, they follow from (19) on elimination of the coefficients aj,6j. Here d'^f is the 
symmetric differential of /. Notice that d'^f and the six quadratic expressions on the 
right hand side of (19) are nothing but second fundamental forms of the submanifold 
C A^ C P^^ defined by (17). Furthermore, the six fundamental forms on the right 
hand side of (19) are the restrictions to of the second fundamental forms of A^ C P^^. 
Thus, (19) says that has no nontrivial second fundamental forms 'of its own', that is, 
all its second fundamental forms can be obtained as restriction of the second fundamental 
forms of A^. This is an obvious necessary condition for a submanifold C A^ C P^^ to 
be a hyperplane section. In the present case, it is also sufficient. Calculating consistency 
conditions for (19) we obtain the equation for and 6j, 

dttQ = a^uj — 2sdui2, dai = aiUJ + sduu, da2 = a2Uj + sdu22, 

(20) 

dbo = Bqu + sdf, dbi = biu + 2sduo2, db2 = b2UJ + 2sduoi, ds = su. 

Here u = a^duxy + aiduyy + a2duxx, and s is yet another auxiliary function. One can 
verify that du = 0. Equations (19) and (20) constitute an involutive differential system 
for / which characterises symplectic Monge- Ampere equations. 

Once we know that our PDE is of symplectic Monge- Ampere type, there are two ways 
to proceed. The first one is to use the fact that linearisable Monge-Ampere equations 
correspond to special hyperplane sections of A^ such that the corresponding hyperplane 
is tangential to A^, that is, belongs to the dual variety [25, 15]. Written in differential 
form, this simple geometric property gives just one extra condition which can be used 
to express s in terms of aj,6j. The resulting formula is quite complicated, reflecting the 
fact the the dual variety of A^ is a quartic hypersurface defined by a rather cumbersome 
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equation. Explicitly, we have: 

■5 = ~{(^0fuxxfuyy + «l/n^^ + (^ifuyy " (^0^' 1 fu ^ X f U^y ~ O,0O-2 f U^y fuy y + ai(J'2{fuxy ^ l U^X l Uyy) — 
(^O^oi fuxtfuyt "I" fuxy) + (^1^1 ifuxtfuxy fuytfuxx) ~^ '^2 ^2 ( /?iy ( /mj:,, fuxtfuyy)~^ 

O'lboifuxt + + a2&o(/4t + — (^ohfuxtfuyy — o.ob2fuytfuxx + 

aihfuxtfuxx + a2bifuytfuyy +bl- blfuyy - blfuxx - bobifuyt - bobifuxt - &1&2A, J/4 det P, 

here the matrix P was defined at the beginning of the proof. A direct calculation shows 
that the remaining coefficients of the Cotton tensor vanish if and only if the above lin- 
earizability condition holds. This, however, is a rather complicated calculation. 

Another way is to use the fact that, in three dimensions, any non-degenerate symplectic 
Monge-Ampere equation is either linearizable, or S'p(6)-equivalent to one of the three 
canonical forms, 

HessM = 1, Hess M = Mil + M22 + W33, Hessw = Mh + M22 — W33, (21) 

see [40, 3]. Note that the first equation governs improper affine hyperspheres, while the 
last two describe special Lagrangian 3-folds. One can verify by a direct calculation that 
formal linearizations of the equations (21) are not conformally fiat for generic solutions. 
Thus, once again conformal fiatness proves to be equivalent to the linearizability. This 
finishes the proof of Theorem 6. 

As in the previous cases, formal linearizations of integrable equations of the form (5) 
give rise to conformal structures satisfying the Einstein- Weyl property. 

Example 1. The equation e"" = e^^"^ + e"^^ appeared in [25] in the classification of 
integrable PDEs of the form F{uxx,Uyy,Utt) = 0. 
Conformal structure: g = e~"^^(ix^ + e'^y^dy'^ — e~""(it^. 
Covector: uj = 2{utttdt + u^xxdx + Uyyydy) - d{utt + m^x + Uyy). 

Example 2. The equation Uu = ^ + ^ri{uxx)ult appeared in [47] in the classification 
of integrable hydrodynamic chains. Here the integrability conditions reduce to the Chazy 
equation for rj, r]"' + 27777" = 8(7/')^. 

Conformal structure: g = Auxtdxdy — (^^rj'u^^ + s^) dy'^ +2sdydt—dt'^ , here s = ^rju^^ — ^. 
Covector: 

^ = (| UtxV + 4:UxyUf^)uttx + (| uf^v"^ + I uf^T]' — ulyUfJ' — I UxyUp}ri)Utxx 

+ (| U^txV V + i U^txV" - \ Uxyi)Uxxx + {UxyU^^ - \ r])Uxxy - 2 Mtey dy 
[(^xj/'^Jx ~l~ '^V)^txx ~l" 3''? '^tx'^xxx '^t,x^xxy 2ltj^ W^jajj dt. 

This structure is Einstein- Weyl if and only if t] solves the Chazy equation. 
The second main result of this Section is as follows: 
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Theorem 7 Equation (5) is integrable by the method of hydrodynamic reductions if and 
only if the corresponding conformal structure g is Einstein- Weyl on any solution, with the 
covector uj given by (8). 

Proof: 

Given an equation of the form (5), the conformal structure of its formal linearization 
is (7 = gijdx^dx^ , where gij is the inverse of Fij. We will seek a covector u = uj^dx^ in the 
form Uk = Tl^' Uxixix' where Tf^ are certain functions of the second order derivatives of u. 
Imposing the Einstein- Weyl equations (7) and using (5) and its differential consequences 
to eliminate all higher order derivatives of u that contain differentiation by t more than 
once (we use the representation (17); maximum fourth order derivatives of u will occur 
in this calculation), we obtain a set of relations, which have to vanish identically in the 
remaining partial derivatives of u. Thus, equating to zero coefficients at the fourth order 
derivatives of u, we obtain the expression (8) for u. 

With this expression, R — Kg is purely quadratic in the third order derivatives u^ixixi- 
Choosing A in such a way that dt^ term disappears, we get 5 ■ 28 = 140 coefficients at 
these quadratic terms which are third order differential polynomials in /. Note that these 
coefficients are linear in the third order derivatives of /. Their vanishing is equivalent 
to some 35 identities constituting an involutive closed system of third order PDEs for /. 
These are precisely the integrability conditions as derived in [25]. This finishes the proof 
of Theorem 7. 

5 Type 4: systems of hydrodynamic type 

In this section we discuss geometric aspects of integrable systems of hydrodynamic type 
(6), 

A{u)u.j: + B{u)uy + C{u)ut = 0, 

where u = {u^^u^Y is a two-component column vector of the dependent variables, and 
A{u), B{u)^C{u) are 2x2 matrices. It will be assumed that there is a matrix in the 
span of A,B,C that is hyperbolic. This class is invariant under arbitrary changes of 
variables u^,u'^, as well as linear transformations of x,y,t, which constitute the natural 
equivalence group of the problem. If C is non-degenerate, the multiplication by 
brings the system into evolutionary form. The integrability of systems of hydrodynamic 
type was investigated in [24]. In was shown in [46] that 'generic' integrable system of the 
form (6) can be parametrized by generalized hypergeometric functions. The linearized 
system has the form 

A{u)vx + B{u)vy + C{u)vt H = 0, 

where dots denote terms which do not contain derivatives of v. The corresponding dis- 
persion relation is given by the formula 

det{X^ A{u) + X^B{u) + X^C{u)) = 0. 

This is a quadratic form in A which can be represented as {X^,X'^,X^)D{u){X^,X'^,X^y 
where D is a 3 x 3 symmetric matrix. It defines the conformal structure g = g^dx'^dx^ 
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where the matrix of gij is the inverse of D. Recall that 
result is as follows: 



(x, y, t). Our first 



Theorem 8 System (6) is linearizahle by a transformation from the equivalence group 
if and only if it is integrable, and the conformal structure g is conformally flat on any 
solution (note that conformal flatness alone is no longer sufficient for the linearizability). 

Proof: 

Without any loss of generality one can assume that system (6) is represented in the form 



u 



u 




p q 



r s 








(22) 



(multiply by C^^ and use a change of variables u^,u^ to make the matrix C~^A diagonal: 
such diagonalization is always possible in the two-component hyperbolic situation). Here 
the matrix elements a,b,p,q,r, s are functions of u^,u^. The corresponding dispersion 
relation takes the form 



(A^a + A'p + X'){XH + + A^) - qr{\^f 



0, 



(23) 



with the associated matrix 



/ ab 



Diu) 



as+hp 



as+bp 
2 



ps — qr 



\ a+b p+s 
\ 2 2 



a+b \ 
9 \ 



1 / 



This gives rise to the conformal structure g = gijdx^dx^ where the matrix of g^j is the 
inverse of D. We require g to be conformally flat for any background solution u. Calculat- 
ing the Cotton tensor (10) and using (22) and its differential consequences to eliminate all 
higher order derivatives of u which contain differentiation by t, we obtain the expression 
which has to vanish identically in the remaining higher order derivatives of u (maximum 
third order derivatives will appear in this calculation). In particular, requiring that co- 
efficients at the remaining third order derivatives of u vanish identically, we obtain that 
a, b,p, s and qr must be constant (thus, all coefficients of the dispersion relation are con- 
stants). In this case the rest of the Cotton tensor vanishes identically. By a transformation 
from the equivalence group any such system can be brought to the form 



where q is still an arbitrary function of u^,v?. For generic q such systems are neither 
linearizable nor integrable. Imposing the integrability conditions as derived in [24] we 
obtain just one extra constraint, (lng)ui„2 = 0, which is equivalent to the existence of 
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a change of variables — j- ip^{u^), — if'^{u^) bringing system (22) to a constant 
coefficient form. This finishes the proof of Theorem 8. 

Formal linearizations of integrable systems of the form (6) give rise to conformal struc- 
tures satisfying the Einstein- Weyl property. 

Example 1. Let us consider the system 

which was obtained in [23] as the ffist order form of the Boyer-Finley equation: indeed, 
the expression p = (m^)^ + {u^Y satisfies the equation 

ptt = Alnp, 

A = d"^ + dy. In this case the dispersion relation is 

P 

The corresponding Einstein- Weyl structure takes the familiar form [55]: 
Conformal structure: g = pi^dx^ + dy"^) — dt^. 

Covector: lo = 2^dt. 

p 

Example 2. Hamiltonian systems of hydrodynamic type can be represented in the form 

Vt + {H,)y = 0, {H^)y = 0, (24) 

here H{v,w) is the Legendre transform of the Hamiltonian density. The dispersion rela- 
tion, 

(A^if^,„ + A^)(A^if^^ + A^) — {X^Hy^Y = O5 
gives rise to the conformal structure 

g = {dy — Hwwdx — H^^dtY — AH^^dxdt. 

The corresponding covector u is given by the formula (8), 

— (^Ay + {Hyy)x + {H^w)i^dy + (^HyyAy + 2Hl^{Hyy)y + Hyy{{Hyy)x + {Hw^)t)^dt, 

here A = H^^Hyj^ — 2H^^. One can verify that g, u satisfy the Einstein- Weyl constraints 
if and only if the potential H{v,w) satisfies the set of integrability conditions as derived 
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in [23] based on the method of hydrodynamic reductions: 



This system is in involution. It was shown in [29] that its 'generic' solution is given by 
the formula 

H{v, w) = Z{v + w) + eZ{v + ew) + e^Z{v + e'^w); (26) 

here e = e^'^*/^ and Z"{s) = (^{s) where ( is the Weierstrass zeta-f unction: (' = —p, {p'Y = 
4:p^ — (equianharmonic case (72 = 0). Degenerations of this solution correspond to 

H{v, w) = -v'^({w), H{v, w) = {v + w) \n{v + w), 
as well as the following polynomial potentials: 

a 1 
H(v,w) = v'^w'^, H(v,w) = vw'^ + —w^, H{v,w) = vw + -w^. 

5 6 

The second main result of this Section is as follows: 

Theorem 9 System (6) with nonconstant dispersion relation is integrable by the method 
of hydrodynamic reductions if and only if the corresponding conformal structure g is 
Einstein-Weyl on any solution, with the covector u given by (8). 

Proof: 

We will give two proofs of this result. The first one is computational, based on the 
explicit calculation of the Einstein-Weyl constraints, and the integrability conditions as 
derived in [24]. The second proof utilises the fact that any two-component integrable 
system of hydrodynamic type possesses a dispersionless Lax pair [24]. We demonstrate 
that the existence of a Lax pair implies that the Weyl connection D, specified by the 
conformal structure g = gijdx''dx^ and the covector (8), possesses a two-parameter family 
of null totally geodesic surfaces, which is known to be equivalent to the Einstein-Weyl 
property [10]. 

The first proof can be summarized as follows. We consider system (6) represented 
in the form (22). The associated conformal structure is (7 = gijdx^dx^ where gij is the 
inverse of D. We will seek a covector u = Ukdx^ in the form Uk = T^jU-'^i where T^j 
are certain functions of u^,u^. Imposing the Einstein-Weyl equations (7) and using (22) 
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and its differential consequences to eliminate all higher order derivatives of u that contain 
differentiation by t (maximum second order derivatives of u will occur in this calculation), 
we obtain a set of relations which have to vanish identically in the remaining partial 
derivatives of u. Thus, equating to zero coefficients at the second order derivatives of m, 
we obtain the expression (8) for uj. Equating to zero the remaining terms we obtain 15 
relations providing all second order partial derivatives (in of the coefficients of the 

dispersion relation, that is, of a, s and qr. These relations constitute part of the total 
set of 16 integrability conditions as derived in [24]. To recover the missing condition we 
calculate the compatibility conditions for the 15 relations at hand. This leads to the two 
cases: 

Case 1: all coefficients of the dispersion relation are constant, in this case all of the 
15 relations are satisfied identically. Any such system can be brought to the form u\ = 
qUy, Uy = qui, see the proof of Theorem 8. Although the Einstein- Weyl property is 
trivially satisfied, the system does not need to be integrable for generic q{v},v?). This is 
why we need to eliminate the case of constant dispersion relation from the statement of 
Theorem 9. 

Case 2: one of the 15 relations, namely, the expression for the mixed derivative {qr)uiy2, 
splits into two separate expressions for qu^u'^ and ru^y2. This gives all of the 16 integrability 
conditions as derived in [24] based on the method of hydrodynamic reductions, thus 
finishing the first proof. 

As in Theorem 5, there exists another demonstration that the integrability implies 
Einstein- Weyl property. It is based on the fact that any integrable system of hydrody- 
namic type, which we again assume represented in the form (22), possesses a dispersionless 
Lax pair 

St = f{Sy, n\ u^), Scc=g{Sy, n\ u^). (27) 

This means that the consistency condition Stx = S^t is equivalent to the system (22). 
Differentiating (27) by y and setting 5*^, = A we obtain 

At = fxXy + fu^ul + /„2My, = gxXy + gu^ul + Qu^uI- (28) 

With this system we associate the pair of vector fields 

which live in the extended four-dimensional space with coordinates x,y,t,X. Note that 
the compatibility condition Xtx = Xxt is equivalent to the commutativity of these vec- 
tor fields. Projecting the two-parameter family of integral surfaces of the distribution 
spanned by X, Y to the space of independent variables x, y, t we obtain a two-parameter 
family of totally geodesic null surfaces of the Weyl connection D. This can be seen as 
follows. Equations (28) mean that the vectors X, Y are tangential to the hypersurface 
A = X{x, y, t). Projecting X and Y to the space of independent variables x, y, t we obtain 
two vector fields 

V- ^ 9 ^_ d d 
ot oy ox oy 
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which commute if and only if A satisfies the equations (28). It remains to show that X and 
Y form a null distribution (that is, tangential to the null cones of the conformal structure 
g), and that the covariant derivatives Dj^X, Vij^Y, DyX, DyF belong to the span of 
X, Y. Equivalently, one can introduce the covector 6 = gxdx + dy + f\dt which annihilates 
X, Y, and verify that 9 is null, and that Bj^O A9 = ByO A9 = 0. This follows from the 
equations satisfied by the functions f{\,v},u^) and g{\,u^,u'^) as derived in [24]: 

a~b ' a-b ' 

^ fful [g^2 (^(^1 )+9^ig^2 {{a-b)p^i +(.s-p)a^i -m^a )+ga„iggi] 

_ 9^2 [g^i (g(a„2 -b„2 )+(fc-a)g„2 )+g„ig„2 {{h-a)s^2 +{v-s)h^2 -qb^i )+rb^2 glj] 

{b-a)qgl^ • 

Note that fx and gx satisfy the relation 

{agx +P + f\){hgx + s + fx) - qr = 0, 

which means that the covector 6 is null in the conformal structure defined by the dispersion 
relation (23). This finishes the proof of Theorem 9. 



6 Further integrable examples in 3D 

In this section we collect miscellaneous examples of dispersionless integrable systems in 
3D related to Einstein- Weyl geometry, which do not fit into the classes discussed above. 

Example 1. The following system was derived by Manakov and Santini [43] as a two- 
component generalisation of the dKP equation: 

Uxt - Uyy + iuUx)x + V^U^y - VyU^x = 0, V^t - Vyy + UV^x + V^V^y - VyV^^ = 0. 

Its formal linearization results upon setting n — n + eui, v ^ v + evi which gives 

Lui H = 0, Lvi^ = 0, 

here L = dxdt—dy + {u—Vy)d'^+Vxdxdy, and dots denote terms which do not contain second 
order derivatives of ui and vi. The symbol of L gives rise to the conformal structure 

g = (dy — VxdtY — 4:{dx — (u — Vy)dt)dt, 

which satisfies the Einstein- Weyl equations with the corresponding covector u given by 
the formula (8): 

W = -Vxxdy + {4:Ux - 2Vxy + VxVxx)dt. 
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This Einstein- Weyl structure was obtained in [17]. 

Example 2. The following system was proposed by Bogdanov [4] as a two-component 
generalisation of the BF equation: 

{e~'^)tt = mt4>xy - mr,(j)yt, mue'^ = m^myt - rritm^y. 

Its formal linearization results upon setting — )■ + e0i, m — m + errii which gives 

L(f)i H = 0, Lmi H = 0, 

here L = e~'^df — m^dydt + rritdxdy, and dots denote terms which do not contain second 
order derivatives of 0i and mi. The symbol of L gives rise to the conformal structure 

g = {m.j.dx + rritdtY + Ae~'^mtdxdy, 

which satisfies the Einstein- Weyl equations with the corresponding covector u given by 
the formula (8): 

UJ = — 5" — z — [nix dx + rrit dt) + 2 — — dy. 
\mf nit J f^t 

Example 3. The equation 

mlmtt = m^niyt - mtnir^y (29) 

was obtained in [4] as a reduction of the two-component BF system from Example 2: 
set = — plnm^. Since this equation fits into the class of quasilinear wave equations 
discussed in Sect. 3, it gives rise to the Einstein- Weyl geometry with the corresponding 
conformal structure 

g = {m^dx + rritdtY + ^m^^^dxdy, 
and the covector u given by (8): 

UJ = (2p -|- 1) — Trlm^ dx + mt dt) + 2 — — dy. 

mf mt 

Note that these g and u result, via the substitution = — plnm^, from the analogous 
formulas of Example 2. It was pointed out in [5] that, applying to this equation the 
contact transformation defined as 

(x, y, t, m, m^, niy, nit) -> {x, y, r, F, F^, Fy, F^) 

where 

1 , ^ m — tnit ^ nix ^ i^v ^ , — if^t 

r=-lnmt, F= F^ = Fy = Fr = -2t^/nit 



Jnit J"<^t yf^t \/^t 



one obtains the equation derived by Dunajski and Tod in the context hyper-Kahler metrics 
with conformal symmetry [20], 

{Fy + Fy,){Fx - Fxr) -{F- F,,)Fxy = Ae^'\ (30) 



23 



It was shown in [20] that this equation gives rise to the Einstein- Weyl structure 



g = (Fdr + F^dx - Fydy - dF^f + We^^^dxdy, 

, , {2 + 4p){F, - F,r)dx + {2 - 4p){Fy + Fyr)dy 
u - Apdr + y-^^ . 

One can verify that the Einstein- Weyl structure of the equation (30) can be obtained 
from that of the equation (29) by applying the above contact transformation (and an 
appropriate rescaling). Note that in the last case the covector u is no longer given 
by our formula (8) (unless p = 0, in which case the equation becomes translationally 
invariant). The reason for this is that, although the Einstein- Weyl property is clearly 
contact-invariant, this is not the case for our formula for u. We recall that, given a 
conformal structure g, the problem of reconstruction of the corresponding covector u 
from the Einstein- Weyl constraints is far from trivial [22]. What simplifies this problem 
in our case is that u should be given by a 'universal' formula depending on finite order jets 
of a solution u. In any case, a fully contact-invariant approach to dispersionless integrable 
systems in 3D is yet to be developed. 

7 Integrability in 4D and self-duality 

There are very few classification results in 4D. Here we consider the case of symplectic 
Monge-Ampere equations, that is, equations represented as linear combinations of mi- 
nors of the Hessian matrix of a function u{x^ , x^ , x^ , x^) , which constitute a subclass of 
equations of type 3 (in 4-dimensional situation). 

Theorem 8 [15] Over the field of complex numbers any integrable non- degenerate sym- 
plectic Monge-Ampere equation is Sp{8)- equivalent to one of the following normal forms: 

1. Mil — U22 — u^i — W44 = (linear wave equation); 

2. ni3 -|- U2A + U11U22 — = (second heavenly equation [49]); 

3. ni3 = ^12^44 — ^14^24 (modified heavenly equation); 
4- 'Wi3'W24 — W14M23 = 1 (first heavenly equation [49]). 

5. Mil + ''^22 + U13U24 — U14M23 = (Husain equation [34]); 

6. aui2U^i + /3iii3^i24 + 7^14^23 = (general heavenly equation), a -|- /3 -|- 7 = 0. 

These heavenly-type equations are known to possess Lax pairs of the form [Xi, X2] = 
where Xi,X2 are parameter-dependent vector fields which commute modulo the corre- 
sponding equations [49, 34, 15, 42]: 

Second heavenly equation 1x13 + ^24 + U11U22 — ^12 = 0: 

Xi = ^4 + ^11^2 - Ui2di + A9i, X2 = ^3 - ^12^2 + U22di - \d2. 
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Modified heavenly equation Ui^ = — U14U24: 



Xi = uud2 - uud^ + Xd] 



X2 



83 + U44d2 - U2AdA + ^84,. 



First heavenly equation M13M24 — M14W23 = 1: 



X2 = -M23C^4 + M24C^3 " Xd2- 



Husain equation Un + ^22 + ^ii3^^24 — ^14^23 = 0: 



Xi = 82 + ui^di - uiidni + Xdi 



X2 = di- U23dn4: + U24d3 - Xd2. 



General heavenly equation aui2U34 + /3mi3M24 + 7W14M23 = 0: 

X, = d,- + 7A(9i - ^^3), X2 = ^84 -82 + pX{d2 





M34 



One can show that the distribution spanned by Xi and X2 is totally null with respect to 
the conformal structure provided by the symbol of formal linearization. The integrability 
of this distribution for any value of the spectral parameter A implies the existence of a 
three-parameter family of null surfaces (a-surfaces) , the property known to be equivalent 
to self-duality Another way to see this is to notice that the above Lax pairs are linear in 



Furthermore, the corresponding formal linearizations can be written in the form 



where the expression XQ — YP is understood as a second order differential operator, 
and dots denote terms containing first order derivatives of v. According to [19, 32] this 
means that the symbols of formal linearizations must be self-dual (in fact, hypercomplex, 
but this property is known to imply self-duality). These and other examples support the 
following conjectures relating the linearizability/integrability of four-dimensional PDEs 
to conformal geometry of symbols of their formal linearizations: 

• A 4D second order dispersionless PDF is linearizable if and only if the corresponding 
conformal structure is conformally fiat on any solution. 

• A 4D second order dispersionless PDF is integrable if and only if the corresponding 
conformal structure is (anti) self-dual on any solution. Since the equations of self- 
duality are known to be integrable by the twistor construction [48], this supports 
the evidence that solutions to integrable PDFs must carry integrable geometry. 

A convenient approach to the integrability of four-dimensional PDFs is based on the 
requirement that all their travelling wave reductions to 3D must be integrable. This 
necessary condition turns out to be very strong indeed, and in many cases is already 
sufficient for integrability [15]. Since the integrability in 3D is related to the Finstein- 
Weyl condition, this gives yet another confirmation of the well-known fact that symmetry 
reductions of the self-duality equations lead to Finstein-Weyl geometry [35, 8, 9]. 



A: 



Xi = X + Ay, X2 = P + XQ. 



[XQ - YP)v + ■ ■ ■ = 
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8 Concluding remarks 



(a) We have characterized the integrabihty of several classes of dispersionless PDEs in 
3D by the Einstein- Weyl property of their formal linearizations. It should be emphasized 
that the four types of equations under consideration, although invariant under certain 
equivalence groups, are clearly not contact/point invariant. On the other hand, our 
characterization is manifestly contact-invariant: the formal linearization of a nonlinear 
PDE plays the role of the 'tangent bundle' to the corresponding solution space. As 
such, its geometric properties (in particular, the Einstein- Weyl property) are invariant 
under arbitrary contact (and more generally Lie-Backlund type) transformations. This 
suggests a contact-invariant approach to the dispersionless integrabihty in 3D unifying all 
known examples. We should however warn the reader that, in general, the Einstein- Weyl 
property alone may not be sufficient for the dispersionless integrabihty: as we saw in the 
proof of Theorem 8, the system u] = qUy, Uy = qu"^ has a conformally flat symbol, but is 
not integrable for generic q{u^,u'^). 

(b) In the case of higher order dispersionless PDEs in 3D, the symbol of formal lineariza- 
tion defines a generalized conformal structure which supplies each solution with a field of 
algebraic null cones. In the spirit of [10], one can define the 'generalized Einstein- Weyl' 
property by the requirement of the existence of a symmetric connection which preserves 
this generalized conformal structure, and possesses a two-parameter family of null totally 
geodesic surfaces. We expect that our results will carry over to this more general context. 

(c) Attempts to extend our results to dimensions higher than four meet an immediate 
difficulty: all known integrable (non-linearizable) examples in dimensions five and higher 
have degenerate symbols. Thus, their solutions do not carry any conventional 'geometry'. 
One possible way to proceed is to require that non-degenerate travelling wave reductions 
of such PDEs to 3D/4D give rise to Einstein- Weyl/self-dual geometries. 

(d) Einstein- Weyl structures are known to be related to third order ODEs satisfying 
the so-called Wiinschmann and Cartan relations [11, 54, 45]. Thus, any solution of a 
dispersionless integrable system has a third order ODE naturally associated with it (to be 
precise, a point equivalence class of ODEs). The (complexified) space of the dependent 
and independent variables of this ODE, which can be identified with the space of null 
totally geodesic surfaces of the Einstein- Weyl structure, is known as the 'minitwistor' 
space [33, 53]. Solutions of this ODE correspond to curves of totally geodesic surfaces 
passing through a fixed point: these form a three-parameter family, and are known as 
'twistor curves'. It would be of interest to develop an efficient procedure to explicitly 
calculate this ODE for multi-phase solutions provided by the method of hydrodynamic 
reductions. 
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